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Preface 


In  this  thesis,  I  have  analyzed  the  effects  of  three  end 
conditions  on  laminated  plates  made  of  orthotropic  plies. 

The  behavior  of  orthotropic  plates  with  a  single  fiber  orien¬ 
tation  were  also  considered  to  give  the  reader  a  point  of 
reference  as  some  work  in  this  area  has  already  been  accom¬ 
plished.  The  point  of  Saint  Venant’s  validity  was  found  in 
each  case  in  order  to  show  the  region  where  this  principle 
can  be  applied. 

Completion  of  this  study  has  required  many  hours  of 
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study.  I  would  also  like  to  acknowledge  Capt.  Jerry  Stinson 
for  his  helpful  comments  and  guidance  in  the  use  of  the  CDC 
Cyber  74  Operating  System.  I  must  further  acknowledge  the 
comments  and  suggestions  made  by  my  thesis  committee  members. 
Dr.  Ernest  Dorko  and  Dr.  Peter  Torvik.  As  in  any  thesis  work, 
the  success  and  often  the  outcome  depend  largely  on  the  guid¬ 
ance  and  assistance  provided  by  one's  thesis  advisor.  A  mere 
acknowledgement  cannot  express  the  debt  of  gratitude  I  owe  my 
thesis  advisor,  Dr.  Anthony  Palazoto  for  his  valuable  advice 
and  direction  given  throughout  this  work.  Finally,  I  must 
acknowledge  my  wife,  Linda,  for  her  encouragement  and  support 
through  the  many  months  of  work  represented  here. 
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Abstract 

The  structural  stiffness  properties  of  composite  lamin¬ 
ated  plates  are  experimental  information  needed  in  order  to 
design  aircraft  components  manufactured  with  composite 
materials.  A  common  test  used  to  obtain  the  necessary  moduli 
is  to  apply  a  tension  load  on  the  ends  of  the  specimen. 

Little  study,  if  any,  has  been  reported  relating  the  applica¬ 
tion  of  the  Saint  Venant ' s  principle  within  a  multi-ply 
laminated  plate. 

Previously,  researchers  have  shown  that  routine  applica¬ 
tion  of  Saint  Venant' s  principle  in  problems  involving  a 
single  composite  lamina  is  not  generally  justified.  This 
study  offers  justification  for  using  Saint  Venant 's  principle 
in  commonly  laid-up  graphite  epoxy  plate  specimens.  The 
boundary  conditions  considered  are  attempts  at  modeling 
varying  end  supports  for  an  in-plane  loading.  The  analysis 
has  been  carried  out  using  a  finite  element  model  considering 
three  aspect  ratios:  3,  5,  and  10. 


INVESTIGATION  OF  SAINT  VENANT * S  PRINCIPLE 


AS  RELATED  TO  LAMINATED 
COMPOSITE  PLATES 

I .  Introduction 

This  thesis  is  an  investigation  of  Saint  Venant's  prin¬ 
ciple  in  laminated  composite  plates  by  use  of  a  finite  element 
model.  Saint  Venant's  principle  has  played  a  fundamental 
role  in  the  application  of  classical  elasticity  theory  to 
problems  of  practical  interest.  The  statement  of  Saint  Venant's 
principle  does  not  lend  itself  to  a  concise  formulation  which 
can  easily  be  translated  into  mathematical  terms.  The  regions 
of  validity  of  this  principle  in  isotropic  materials  has  been 
extensively  studied  and  is  readily  available  in  literature 
[1] .  However,  the  conquest  of  space  and  the  rise  in  energy 
costs  have  made  the  Air  Force,  as  well  as  the  aerospace  in¬ 
dustry,  very  conscious  of  overall  vehicle  weights.  This  has 
lead  to  the  use  of  lightweight  laminated  composite  materials 
as  primary  load  bearing  structures.  Advanced  composite 
material  use  has  prompted  numerous  studies  of  the  behavior 
of  orthotropic  and  anisotropic  materials  using  either  classi¬ 
cal  elasticity  theory  or  experimental  procedures  to  determine 
material  properties.  The  use  of  the  Saint  Venant  assumption 
is  basic  to  the  measurement  of  these  material  properties  as 
well  as  to  structural  design. 
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Problem  Statement 

This  investigation  used  the  finite  element  analysis  and 
classical  laminated  plate  theory  to  model  the  stress  gradients 
throughout  plates  with  aspect  ratios  (i.e.  length  to  width) 
of  three,  five,  and  ten.  This  modeling  technique  allows  the 
monitoring  of  stress  gradient  propagation  at  discreet  points 
within  the  plate.  As  the  decay  of  the  stress  gradient  would 
logically  depend  on  the  loading  conditions,  different  loading 
conditions  are  considered.  Detailed  discussion  of  the  loading 
conditions  is  found  in  the  numerical  analysis  section.  In 
addition,  the  material  properties  have  a  great  bearing  on  the 
decay  lengths.  The  anisotropic  behavior  of  laminated  plates 
is  directly  related  to  the  ply  orientations.  In  order  to 
ascertain  the  effects  of  ply  orientation  (0,  ±45,  90)s, 

(0,  ±45)s,  (±45)s,  and  (0,  90)s  laminated  plates  are  analyzed. 
The  individual  lamina  characteristics  were  also  investigated 
for  Odeg,  90deg,  and  45deg  orientations.  The  analysis 
utilized  the  Control  Data  Corporation  6600/Cyber  74  computer 
system. 

Throughout  this  study  the  material  properties  for  graphite 
epoxy  listed  by  Ashton,  e_t  al_.  [2]  were  used  in  calculating 
stresses.  These  properties  are  shown  in  Table  I. 

Background 

An  extended  discussion  of  some  of  the  theoretical  formu¬ 
lations  and  further  definitions  of  the  Saint  Venant's  principle 
is  presented  in  a  subsequent  section.  For  purposes  of  the 
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TABLE  I 

Graphite  Epoxy  Properties 


Elastic  Constants 

*  Allowable  Stresses 

t 

.  1  in 

axt 

70 

ksi 

E11 

18.5  x  106 

psi 

°xc 

110 

ksi 

E22 

1.54  x  106 

psi 

V 

6.5 

ksi 

G12 

.85  x  106 

psi 

ayc 

15 

ksi 

v12 

.25 

axy 

8.5 

ksi 

|NOTE:  *Used  as  input  parameter  for  OPSTAT  only. 


present  discussion,  Saint  Venant ' s  principle  is  assumed  to  be 
the  ability  to  replace  complicated  end  loads  with  statically 
equivalent  loads.  The  errors  caused  by  this  assumption  will 
die  out  or  decay  at  a  distance  sufficiently  far  from  the  ends. 
This  definition  closely  parallels  those  given  in  elementary 
solid  mechanics  texts  such  as  reference  [3] .  The  interpre- 
tation  of  this  principle,  as  well  as  its  application,  is  not 
entirely  clear.  The  justifications  for  the  use  of  this 
assumption  are  largely  empirical;  however,  a  large  class  of 
problems  are  idealized  mathematically  using  this  assumption. 
Many  of  these  problems  would  be  difficult,  if  not  impossible, 
to  solve  using  a  classical  elasticity  approach.  An  example 
of  a  difficult  problem  which  has  been  simplified  by  the  use 
of  this  principle  is  the  simple  tension  test  in  which  the  ends 
of  the  specimen  are  clamped.  The  Saint  Venant 's  principle  is 
invoked  by  assuming  a  uniform  force  distribution  sufficiently 
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far  from  the  clamping  devices. 

It  is  evident  from  the  many  formulations  and  examples 
given  in  the  literature  that  the  decay  length  or  distance 
sufficiently  far  from  the  end  forces  in  a  body  depends  on  the 
boundary  conditions,  body  geometry,  and  material  properties. 

It  is  reasonable  to  assume  that  the  region  of  validity  for 
the  Saint  Venant  assumption  will  be  different  for  isotropic, 
orthotropic  and  anisotropic  materials. 

Using  classical  elasticity  theory,  Horgan  [4]  has  investi¬ 
gated  the  decay  of  end  effects  in  anisotropic  materials  and 
has  concluded  that  the  end  effects  decay  slowly  in  highly 
anisotropic  materials.  The  slow  decay  of  stresses  with 
distance  from  the  boundary  of  a  body  under  the  influence  of 
self-equilibrating  forces  gives  rise  to  questions  about  the 
measurement  of  elastic  material  properties.  Folkes  and 
Arridge  [5]  have  shown  experimentally  that  end  effects  in 
highly  anisotropic  polymers  are  very  slow  to  decay  and  can 
give  rise  to  erroneous  material  property  measurements. 

Choi  and  Horgan  [6],  again  using  classical  theoretical 
elasticity,  have  investigated  the  decay  lengths  in  an  ortho¬ 
tropic  material.  Further  detail  of  the  mathematical  solution 
employed  by  Choi  and  Horgan  can  be  found  in  Appendix  D.  This 
analysis  used  a  series  representation  of  the  solution  to  show 
that  the  decay  lengths  along  the  fiber  axis  of  an  orthotropic 
composite  material  is  considerable  in  small  aspect  ratio  plates. 

The  question  of  applications  to  multilayered  composite 
plates  of  similar  aspect  ratios  is  prompted  and  is  the  subject 
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of  this  investigation.  Also,  questions  are  raised  concerning 
orientations  other  than  those  along  the  fiber  axis. 

The  study  of  off-axis  behavior  of  composite  materials  has 
been  reported  by  Pagano  and  Halpin  [7].  This  study  indicates 
that  significant  shear  and  bending  effects  are  present  in 
tension  tests  which  incorporate  clamped  ends.  The  experimen¬ 
tal  results  show  large,  non-uniform  displacements  near  the 
clamped  edges  and  more  uniform  displacements  near  the  center 
of  an  off-axis  plate.  Pagano  and  Halpin  have  shown  large 
shear  coupling  effects  in  the  off-axis  specimens.  Jones  [8] 
has  commented  that  these  results  indicate  that  a  much  longer 
gage  length  is  needed  in  the  off-axis  tension  specimen  in 
order  to  reach  a  region  in  which  Saint  Venant's  principle  may 
be  applied.  Pagano  and  Halpin  have  further  stated  that  a 
possible  alternative  is  to  test  bidirectional  (±e)s  specimens 
for  material  properties. 

Whitney  [9]  has  shown  that  material  properties  measured 
in  tensile  tests  of  symmetric  orthotropic  laminates  may  be  a 
function  of  the  stacking  sequence.  He  further  notes  that  if 
end  conditions  due  to  clamping  are  neglected  in  off-axis 
tensile  tests,  erroneous  values  for  Exx  will  result.  The 
magnitude  of  the  error  is  inversely  proportional  to  the  aspect 
ratio  of  the  plate.  Wu  and  Thomas  [10]  have  further  pointed 
out  that  the  clamped  end  condition  induces  a  moment  which  is 
not  a  part  of  the  Saint  Venant  effect.  This  further  compli¬ 
cates  the  experimental  study  of  this  problem.  Rizzo  [11]  , 
using  a  finite  element  model,  has  pointed  out  that  the  rotation 


of  the  specimen  near  the  clamp  is  a  function  of  the  specimen 
aspect  ratio  and  the  rigidity  of  the  clamp  as  well  as  material 
properties.  He  has  concluded  that  rotating  clamps  would  help 
reduce  some  of  the  end  effects,  and  aspect  ratios  greater 
than  six  are  needed.  Richards,  et  al.  [12 j  have  confirmed  the 
findings  of  Rizzo  using  another  finite  element  model.  This 
model  utilized  two-dimensional,  constant  strain  triangles 
similar  to  the  one  used  in  this  study. 

Although  much  can  be  found  in  the  literature  about  the 
orthotropic  composite  plate  under  end  loads,  little  has  been 
said  about  the  laminated  plate  under  similar  loads.  It  is 
hoped  that  this  study  will  help  bridge  that  gap. 
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1 1 .  Theory 


The  theoretical  considerations  in  this  study  may  be  summed 
up  by  answering  the  following  questions.  What  is  the  general¬ 
ized  Saint  Venant* s  principle?  Where  is  this  principle  valid? 
Does  the  model  used  correctly  transfer  the  load?  When  does 
laminated  plate  theory  apply?  The  first  two  questions  are 
discussed  in  the  following  subsection.  The  last  two  questions 
have  to  do  with  the  implications  of  a  two-dimensional  model 
and  are  discussed  in  a  separate  subsection. 

Saint  Venant’s  Principle 

In  1855  Saint.  Venant  stated  the  principle  of  elastic 
equivalence  of  statically  equal  systems  of  forces.  This  prin¬ 
ciple  was  stated  in  justification  of  a  result  which  applied 
rigorously  to  a  load  condition  other  than  the  one  under  study. 
The  original  statement  of  the  principle  has  been  greatly 
generalized  to  include  all  types  of  loading  systems.  Because 
of  the  fundamental  nature  of  this  principle  in  solid  mechanics, 
further  study  of  its  application  in  orthotropic  and  anisotropic 
materials  is  needed,  but  a  solid  understanding  of  the  basic 
concepts  related  to  the  principle  is  required. 

One  formulation  of  Saint  Venant 's  principle,  attributed  to 
Zanaboni  by  Fung  [13]  >  considers  strain  energy  density  as  a 
measure  of  the  region  of  validity  in  an  arbitrary  body.  This 
density  is  shown  to  decrease  in  cross  sections  farther  removed 
from  the  surface  where  a  generalized  force  distribution  anomoly 
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occurs . 


If  one  is  sufficiently  removed  from  the  force  ano- 
moly,  the  strain  energy  density  is  negligible.  The  region  of 
validity  of  Saint  Venant's  principle  under  this  formulation 
is  that  region  where  the  strain  energy  density  to  the  force 
anomoly  is  negligible.  The  Zanaboni  formulation  is  now  as 
follows : 

Let  Sj  and  S2  be  two  nonintersecting  sections  of  a  body 
outside  the  sphere  shown  in  Fig  1.  If  the  section  S2  lies 
at  a  greater  distance  than  the  section  from  the  sphere  in 
which  a  system  of  self-equilibrating  forces  acts  on  the  body, 
then  the  strain  energy  density  caused  by  surface  tractions 
at  R2  is  less  than  the  strain  energy  density  caused  by  sur¬ 
face  tractions  at  R^ ,  where  R1  and  R2  are  the  surface  areas 
at  cross  sections  and  S2,  respectively. 

The  significance  of  this  result  is  that  at  some  point 
sufficiently  far  from  the  sphere,  the  effects  of  the  force 
system  is  negligible.  The  effect  of  the  true  end  condition 
load  is  unimportant  if  Saint  Venant's  principle  is  applicable. 
This  can  be  further  generalized  to  justify  assuming  a  uniform 
edge  loading  where  a  non-uniform  edge  loading  occurs. 

Another  formulation  of  Saint  Venant's  principle  due  to 
Von-Mises  is  presented  to  give  further  insight;  details  can 
be  found  in  Fung.  Von-Mises  considered  four  examples  of 
forces  on  the  boundary  of  a  half  space.  The  forces,  F, 
acted  on  a  small  circle  of  diameter,  e,  as  shown  in  Fig  2. 

The  order  of  magnitude  of  the  largest  stress  component  at 
an  arbitrary  point  which  lies  a  distance  R  from  the  circle 
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is  =  F/R  for  case  a»  and  on  the  order  of  (e/R)oQ  for 

cases  b  through  d  where  aQ  is  the  uniform  stress  value  that 
would  be  assumed  at  a  suitable  distance.  This  leads  to  the 
conclusion  that  end  effects  and  surface  stress  variations 
become  negligible  at  some  point  sufficiently  far  from  the  edge. 

Another  example  of  the  decrease  in  the  end  condition  force 
has  been  demonstrated  by  Timoshenko  and  Goodier  [1]  for  an 
isotropic  plate  using  a  classical  elasticity  series  solution. 
This  example  shows  a  plate  under  a  single  point  compressive 
load  P  at  the  center  of  the  short  ends.  The  stress  at  cross- 
sections  along  the  plate  length  are  shown  to  uniformly 
approach  P/A,  where  A  is  the  cross-sectional  area  at  a  dis¬ 
tance  equivalent  to  the  width  of  the  plate.  At  this  distance, 
the  load  could  be  assumed  equal  to  a  uniform  stress,  oq  =  P/A. 
This  example  is  considered  as  justification  in  using  a  simple 
equivalent  stress  for  other  loading  conditions  by  invoking 
Saint  Venant's  principle  at  a  particular  distance  [4].  In 
general,  this  is  not  true  for  orthotropic  composite  materials, 
as  shown  by  Choi  and  Horgan.  Thus,  the  question  is  raised 
pertaining  to  the  region  of  validity  in  laminated  composite 
plates.  The  study  of  this  region  has  obvious  utility  in 
design  and  measurement  of  material  properties. 

The  author  has  found  several  practical  descriptions  of 
Saint  Venant's  principle  in  the  literature  of  which  two  are 
mentioned  here.  The  first  is  due  to  Filonenko-Boradich  [14] 
and  is  repeated  here:  "At  points  in  a  solid  which  are  suffi¬ 
ciently  far  from  the  surfaces  of  application  of  external  loads, 
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stresses  depend  very  slightly  on  the  particular  manner  in 
which  the  loads  are  applied.”  One  drawback  in  this  defini¬ 
tion  of  Saint  Venant’s  principle  is  the  lack  of  load  defini¬ 
tion.  A  better  definition,  and  the  one  which  has  been  used 
in  this  study,  has  been  offered  by  Sokolnikoff  [15] .  His 
statement  of  Saint  Venant's  principle  is  as  follows: 

If  some  distribution  of  forces  acting  on  a  portion 
of  the  surface  of  a  body  is  replaced  by  a  different 
distribution  of  forces  acting  on  the  same  portion  of 
the  body,  then  the  effects  of  the  two  different  dis¬ 
tributions  on  the  parts  of  the  body  sufficiently  far 
removed  from  the  region  of  application  of  the  forces 
are  essentially  the  same,  provided  that  the  two  dis¬ 
tributions  of  forces  are  statically  equivalent. 

Here  static  equivalence  refers  to  equivalent  resultant  forces 

and  moment.  This  definition  gives  the  necessary  equivalence 

of  force  guidelines. 

The  previous  definitions  and  concepts  have  answered  the 
first  two  questions  raised  at  the  opening  of  this  section. 

The  previous  discussion  has  given  some  insight  into  how  the 
Saint  Venant  assumption  was  originally  employed  and  how  it  has 
been  generalized  into  a  fundamental  building  block  of  solid 
mechanics.  This  study  has  applied  the  same  definitions  to  the 
laminated  composite  plates  to  give  examples  and  to  determine 
the  appropriate  decay  lengths  of  specific  end  condition 
effects . 


Two-Dimensional  Model  Implications. 

The  application  of  Saint  Venant’s  principle  in  an  elastic 
solid  has  to  be  three-dimensional  for  complete  validity.  In 
order  for  a  two-dimensional  model  to  be  valid,  the  state  of 


stress  in  the  third  dimension  must  be  shown  to  have  negligible 
affect  on  the  stress  gradients  in  the  plane  of  the  two  dimen¬ 
sions  modeled.  Filonenko-Boradich  has  stated  that  the  surface 
of  force  application  must  be  small  compared  to  the  overall 
dimensions  of  the  body,  and  the  surface  of  force  application 
must  be  at  least  of  the  same  order  as  the  smallest  body  dimen¬ 
sion  for  validity  of  the  Saint  Venant  assumption.  This  implies 
that  Saint  Venant' s  principle  cannot  be  applied  across  the 
thickness  of  a  plate.  However,  classical  elasticity  theory 
has  modeled  thin  plates  as  a  two  dimension  body  with  great 
success.  The  Timoshenko  and  Goodier  plate  example  previously 
discussed  uses  this  two-dimensional  model  with  very  good  re¬ 
sults.  This  raises  the  question  of  validity  in  the  laminated 
composite  plate  and  plane  stress  modeling  in  such  plates. 

In  the  work  of  Pipes  and  Pagano  [16],  the  interlamina 
stresses  were  studies.  These  stresses  were  found  to  be  signi¬ 
ficant  only  in  regions  near  the  free  edge  approximately  equal 
to  the  laminate  thickness  where  the  thickness  was  at  least  an 
order  of  magnitude  smaller  than  the  other  dimensions.  Justi¬ 
fication  for  the  plane  stress  model  and  for  the  invocation  of 
Saint  Venant *s  principle  in  the  third  dimension  exists  if  that 
dimension  is  small  compared  to  the  other  two  dimensions. 

If  one  examines  classical  laminated  plate  theory  [3],  it 
is  found  that  the  stress  gradient  across  the  lamina  thickness 
is  averaged  out  of  the  problem  and  only  the  two-dimensional 
portion  of  the  problem  is  considered.  If  the  interlamina 
stresses  are  not  significant,  one  might  also  ask  if  the 
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averaging  of  the  stress  across  a  lamina  as  in  classical 
laminated  plate  theory  has  an  appreciable  affect  on  the  vali¬ 
dity  of  solutions.  The  analysis  of  Rizzo  and  that  of  Richards, 
et  al . ,  are  essentially  dealing  with  a  plate  of  a  single  off- 
axis  lamina.  The  two-dimensional  models  used  in  these  studies 
showed  a  high  correlation  with  experimental  results. 

Probably  the  most  attractive  reason  for  using  a  two- 
dimensional  model  is  the  level  of  computational  difficulty 
compared  to  the  three-dimensional  problem.  The  two-dimensional 
solution  has  given  good  results  in  the  isotropic  case  and  in 
the  laminated  composite  examples  cited,  with  less  computation. 
The  level  of  computational  difficulty  in  many  cases  can  be 
directly  tied  to  the  cost  of  computer  time  and  ultimately,  in 
some  cases,  to  the  ability  to  obtain  a  solution  for  a  parti¬ 
cular  problem. 
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III.  Numerical  Analysis 

The  finite  element  method  was  used  to  generate  the  data 
obtained  in  this  thesis.  The  finite  element  models  incorpor¬ 
ated  in  the  numerical  analysis  are  discussed  herein. 

General  Comments 

In  the  previous  sections,  it  was  pointed  out  that  a  point 
exists  sufficiently  far  from  the  end  conditions  where  the 
stress  no  longer  becomes  affected  by  the  loading  conditions. 

In  order  to  locate  this  point,  an  exact  solution  becomes 
necessary.  Since  exact  solutions  are  difficult,  if  not  im¬ 
possible  to  obtain  for  composite  material  problems,  it  is 
useful  to  carry  out  approximate  numerical  solutions  and  place 
a  level  of  confidence  on  the  results  obtaineu. 

The  level  of  confidence  placed  on  a  given  solution  obtained 
due  to  a  particular  finite  element  model  may  be  associated  with 
the  error  involved  in  the  technique  used  to  obtain  the  solution. 
Durocher  and  Palazotto  [17]  have  pointed  out  that  the  sources 
of  error  can  be  classed  as  idealization,  discretization  and 
manipulative  errors.  Replacing  an  anisotropic  or  orthotropic 
material  with  a  homogeneous  continuum  is  a  typical  example  of 
the  idealization  error.  Errors  of  this  type  are  common  to  all 
engineering  models.  The  discretization  error  is  due  to  re¬ 
placing  the  continuous  medium  with  finite  elements.  The 
manipulative  errors  are  associated  with  truncation  and  round 
off  produced  by  computer  programming.  The  latter  sources  of 
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error  are  not  inherent  to  the  engineering  model.  Thus,  a 
convergence  study  is  required  to  guarantee  the  numerical 
solution  approaches  the  exact  or  idealized  solution.  The  con¬ 
vergence  of  the  chosen  numerical  models  will  be  discussed  in 
subsequent  sections. 

The  finite  element  program  used  in  this  analysis  was 
developed  by  Venkayya  [is]  of  the  Air  Force  Flight  Dynamics 
Laboratory.  Classical  laminated  plate  theory  and  the  displace¬ 
ment  method  of  finite  element  is  used  in  this  program.  The 
program  can  be  utilized  to  analyze  plates  with  several  combina¬ 
tions  of  ply  orientations  as  well  as  orthotropic  plates.  This 
study  employed  (0,  ±45,  90)s,  (0,  90)s,  (0,  ±45)s,  and  (±45)s 
oriented  plates  along  with  plates  of  Odeg,  90deg,  and  45deg 
fiber  orientations.  The  program  incorporated  several  standard 
elements,  including  the  constant  strain  triangle  which  was 
exercised  in  this  study.  Relative  percentage  of  ply  orienta¬ 
tions  contained  in  the  laminate  and  element  material  properties 
were  applied  in  calculating  the  elemental  stiffness  matrix. 

The  necessary  characteristics  of  proper  convergence  is 
noted  by  Durocher  and  Palazotto.  Mesh  refinement  has  been 
carried  out  to  enforce  convergence.  The  specific  end  conditions 
will  now  be  addressed  using  Fig  3. 

Parabolic  Tension  End  Condition 

The  parabolic  tension  end  condition  is  an  attempt  to  model 
a  non-uniform  tensile  end  force  in  which  the  edges  of  the  plate 
are  allowed  to  slip.  The  plate  is  loaded  at  the  ends  with  a 
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where  N  is  the  end  force,  N  is  the  normalization  force,  and 
y  is  the  coordinate  across  the  width  of  the  plate  with  zero 
at  the  center.  The  load  is  symmetric  about  the  x-axis  and  the 
y-axis,  thus  only  a  quarter  of  the  plate  need  be  modeled.  The 
material  properties  for  all  the  laminated  plate  orientations 
and  the  Odeg  and  90deg  fiber  oriented  plates  are  symmetric 
about  these  axes  also.  The  exception  to  this  is  the  45deg 
fiber  oriented  plate  which  will  be  discussed  separately.  The 
first  quadrant  meshes  used  for  these  cases  are  shown  in  Fig  4. 

In  order  to  model  the  high  stress  gradients  with  constant 
strain  triangles,  a  more  refined  mesh  is  needed  over  those 
areas.  This  dictates  a  larger  number  of  elements  near  the 
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Fig  4.  Quarter  Plate  Convergence  Study  Meshes  With  Symmetry 
Conditions  for  Parabolic  Tension  End  Conditions 


loads.  The  symmetry  conditions  require  restraints  as  shown 
in  Fig  4  along  the  quadrant  boundaries. 

As  constant  strain  triangles  were  used,  the  centroid 
coordinates  of  the  triangle  are  employed  for  stress  comparison 
and  the  loci  of  points  must  be  established  in  the  analysis  for 
convergence  to  a  particular  solution.  The  end  condition  under 
study  was  also  investigated  by  Choi  and  Horgan  in  their  analy¬ 
sis  of  an  orthotropic  graphite  epoxy  plate.  The  series 
solution  obtained  by  Choi  and  Horgan  was  compared  to  the 
results  in  the  convergence  study. 

Three  mesh  arrangements  were  studied  in  order  to  estimate 
convergence.  Stresses  were  checked  along  the  width  of  the 
plate  at  various  stations  on  the  x  (length)  axis.  The  ox 
stress  across  width  sections  must  average  to  2/3  of  the  stress 
normalization  factor  aQ,  and  oy  stresses  must  average  to  zero. 
This  was  the  case  in  the  three  mesh  arrangements  to  within 
five  percent.  The  fine  mesh  arrangement  shown  in  Fig  4-C 
had  average  stresses  ox  values  within  one  percent  of  2/3  0Q. 

All  three  solutions  compared  favorably  with  the  Choi  and  Horgan 
solutions . 

As  previously  touched  on,  the  45deg  oriented  fiber  plate 
could  not  be  considered  symmetrical  because  of  material  proper¬ 
ties.  This  is  in  contrast  to  the  ±45deg  combination  as  in 
these  laminates  the  plies  act  together  as  one  material.  Figure 
5  depicts  the  boundary  conditions  utilized  for  the  45deg  fiber 
oriented  plate  and  the  mesh  arrangements  which  were  used  for 
the  full  plate  convergence  study.  Two  mesh  arrangements  were 
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Fig  5.  Full  Plate  Convergence  Study  With 
Parabolic  Tension  End  Condition 
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employed  in  this  convergence  study  and  the  results  were 
similar  to  those  obtained  with  the  quarter  plate  mesh  arrange¬ 
ments.  In  both  the  quarter  plate  and  full  plate  mesh  arrange¬ 
ments,  the  loads  were  applied  to  the  end  along  X  =  L. 


Parabolic  Bending  End  Condition 

The  parabolic  bending  end  condition  was  considered  to  find 
the  effects  of  end  rotations  on  the  decay  lengths  for  Saint 
Venant's  assumption.  This  represents  a  pure  bending  load 
which  is  not  found  in  practical  application,  but  this  condition 
gives  insight  to  the  in-plane  rotational  problem  found  in 
rotating  clamping  devices.  The  force  distribution  per  unit 
width  is  given  by 

Nx  =  NQ(5y3/3)  (III-2) 

It  should  be  noted  that  the  term  parabolic  bending  refers  to 
the  shape  of  the  loading  curve  as  shown  in  Fig  6-a.  Again, 

Nv  is  the  end  force,  N  is  the  normalization  force,  and  y  is 
the  coordinate  across  the  plate.  As  the  full  plate  is  modeled, 
all  laminates  and  single  fiber  oriented  plates  were  analyzed 
using  this  mesh  arrangement.  This  mesh  is  the  same  as  that 
applied  to  the  parabolic  tension  case.  Timoshenko  and  Goodier 
have  shown  that  this  loading  condition  produces  a  straight 
line  function  when  the  point  of  Saint  Venant's  principle 
validity  is  reached. 

The  average  ox  stress  for  this  case  is  zero.  All  loci  of 
centroid  points  across  the  width  of  the  plate  averaged  to  zero 


Fig  6.  Mesh  Arrangements  For  Parabolic  Bending 
and  Clamped  End  Conditions 


within  one  percent.  This  load  was  positioned  on  both  ends  of 
the  plate. 

Clamped  End  Condition 

The  final  end  condition  considered  was  the  clamped  end 
condition.  This  is  a  condition  normally  associated  with 
experimental  tension  specimens.  The  shear-moment  coupling  is 
inherent  in  off-axis  plate  tests  and  would  only  be  present  in 
symmetric  laminated  composite  plates  if  rotating  clamping 
devices  were  used  with  an  eccentric  load.  This  loading  was 
investigated  primarily  to  show  the  effects  of  shear  on  the 
finite  element  model  and  to  specifically  point  out  difficul¬ 
ties  encountered  in  using  the  constant  strain  triangle.  Wu 
and  Thomas  call  attention  to  the  fact  that  bending  and  shear 
coupling  effects  are  not  Saint  Venant  type  problems  but  do 
exist  in  real  clamped  end  tension  test.  The  end  condition 
load  functions  are 


Nx  =  V10  '  2XY)  (HI-3) 

and  Nxy  =  Nq(-1  +  Y2)  (1 1 1  -4 ) 

Here  again,  N  is  the  end  force  in  the  x  direction,  N  is  the 

A  v 

normalization  factor  and  X  and  Y  are  the  plate  coordinates. 

N  is  the  end  force  in  the  y  direction  and  provides  the  shear 
xy 

to  compensate  for  the  moment  produced  by  the  2XY  term  in  Nx. 

The  same  full  plate  mesh  was  utilized  in  this  analysis  as  was 
employed  in  the  previous  end  condition  cases. 

The  load  functions  given  above  are  more  simple  than  those 
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used  by  Rizzo  and  can  be  derived  by  an  Airy  type  function 

<j>  =  bxy  +  jy2  +  £xy3  (III-5) 

where  b,  c,  and  d  are  arbitrary  constants.  The  resulting  load 
functions  are  a  uniform  tension  added  to  a  uniform  in-plane 
bending  moment  and  a  parabolic  shearing  force  which  is  re¬ 
quired  for  moment  equilibrium.  These  loads  will  give  insight 
into  the  actions  of  various  plies  and  the  dependence  on  ply 
lay-up  for  strength  measurements  in  laminated  plates. 

Although  the  uniform  loading  in  the  x  direction  and  the 
parabolic  shear  loading  representation  of  the  clamped  end 
problem  is  simplistic,  the  insight  gained  from  this  end  con¬ 
dition  will  give  the  reader  some  feeling  for  the  difference 
between  the  Saint  Venant  decay  length  problem  and  the  moment- 
shear  coupling  problem  in  a  true  clamped  end  plate. 


23 


IV.  Results 


The  results  of  the  finite  element  analysis  of  the  three 
end  conditions  will  be  discussed  in  this  section.  Graphical 
representation  of  the  normalized  stress  dissipation  across  the 
width  of  the  plate  at  specific  stations  along  the  length  of 
the  plate  are  included  in  Appendices  A  through  C.  The  normal¬ 
ized  stress  values  are  given  by 

SIGMAX  =  ax/oo 

(IV-1) 

and  SIGMAXY  =  oxy/oQ 

where  SIGMAX  and  SIGMAXY  are  normalized  stresses,  a  and  a 

a  Ay 

are  the  stress  in  the  x  direction  and  shearing  stress  respec¬ 
tively,  and  aQ  is  the  stress  normalization  factor  computed  by 
dividing  the  resultant  force  normalization  factor,  NQ ,  by  the 
width  of  the  plate. 

The  point  of  Saint  Venant's  principle  validity  is  described 
by  a  constant  uniform  state  of  stress  in  each  case.  The  decay 
length  of  the  a  stresses  were  chosen  by  comparing  the  o 
values  along  a  particular  loci  of  mesh  element  centroids,  across 
the  width  of  the  plate,  with  the  constant  state  of  stress  for 
the  particular  end  condition.  The  station  or  x-coordinate 
which  defined  the  region  of  Saint  Venant's  principle  validity 
was  chosen  as  that  station  where  the  locus  of  stress  values  was 
within  ten  percent  of  the  constant  state  of  stress.  The  ten 
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percent  margin  is  consistent  with  the  mesh  size,  the  program 
round-off  error,  and  the  ability  to  practically  measure  stress 
in  an  actual  application. 

In  all  cases,  the  o  and  a  direct  stresses  acted  similarly 

x  y 

in  that  the  constant  stress  value  was  reached  at  the  same 
station.  The  ox  stresses  were  highest  in  the  Odeg  plies, 
followed  by  the  ±45deg  plies.  The  lowest  ox  stress  values 
were  evidenced  in  the  90deg  ply.  This  is  due  to  the  relative 
stiffness  properties.  Discussion  of  the  individual  end  con¬ 
ditions  follows. 

Parabolic  Tension  End  Condition 

The  Odeg  fiber  oriented  plate  results,  shown  in  Fig  A-l, 
agree  closely  with  the  results  obtained  by  Choi  and  Horgan. 

It  was  found  that  this  case  produced  one  of  the  longest  decay 
lengths  and,  thus,  one  of  the  smallest  areas  for  use  of  the 
Saint  Venant's  principle.  Larger  areas  of  Saint  Venant's 
principle  validity  were  found  for  the  90deg  and  45deg  fiber 
oriented  plates. 

The  laminated  plates  were  investigated  for  three  aspect 

ratios  (length  to  width,  L/W) .  The  constant  uniform  state  of 

stress  for  this  end  condition  averaged  across  the  plies  of  the 
2 

laminate  is  yOQ.  However,  a  different  value  for  the  constant 
uniform  stress  was  found  for  each  ply  according  to  the  relative 
stiffness  of  the  ply.  All  laminas  had  the  same  percent  thick¬ 
ness  for  a  given  ply  lay-up,  and  all  plates  had  the  same  total 
thickness.  Thus,  plates  with  less  laminas  had  thicker  laminas 
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than  plates  with  more  laminas  and  exhibited  lower  stress 
gradients . 

As  shown  in  Table  II,  the  percent  of  the  plate  length 
where  the  stresses  reach  a  value  within  ten  percent  of  the 
uniform  value  is  on  the  order  of  one  plate  width  for  aspect 
ratios  of  five  and  ten,  with  the  exception  of  the  (0,  90)s 
plate.  In  this  plate  the  same  regions  of  constant  uniform 
stress  were  not  exhibited  for  both  plies  in  the  three  and  five 
aspect  ratio  cases,  and  the  decay  length  is  much  longer  in  the 
L/W  =  10  case  as  compared  to  the  other  laminates.  This  is 
attributed  to  the  greater  stiffness  of  the  Odeg  ply  in  tension. 
The  stress  values  differ  in  this  case  by  a  factor  equal  to  the 
ratio  of  E.^  to  E ^ >  which  would  be  expected  for  a  pure  tension 
loading . 

The  aspect  ratio  of  three  had  regions  of  stress  values 
within  ten  percent  of  the  constant  uniform  stresses  for  the 
(±45)s  and  the  (0±45)s  plates.  The  (0,  +45,  90)s  case  exhibited 
dissimilar  ply  actions.  Again,  this  is  attributed  to  the 
stiffness  properties.  The  (0,  90)s  plate  did  not  reach  a  uni¬ 
form  stress  value  for  the  aspect  ratio  of  three.  It  must  be 
pointed  out  that  if  a  margin  of  less  than  ten  percent  were  used 
for  the  aspect  ratio  of  three,  none  of  the  plates  would  reach 
the  uniform  stress  value  for  this  end  condition. 

As  shown  in  Fig  A-3,  the  stress  values  for  the  45deg  fiber 
oriented  plate  had  a  certain  amount  of  variation  about  the 
mean  value.  This  phenomena  is  specifically  pointed  out  since 
all  off-axis  plies  exhibited  this  variation  in  stress  values. 


The  variation  appears  to  be  related  to  the  shear  stresses 
in  the  problem  as  is  later  evidenced  in  the  clamped  end  con¬ 
dition. 


TABLE  II 

Parabolic  Tension  End  Condition 
Uniform  Stress  Station 


Orientation 

L/W  =  5 

L/W  =  10 

0 

i 

.43L 

90 

.  86L 

-- 

45 

.  7  3L 

-- 

(±45) s 

.  63L 

.77  L 

.  87L 

(0,±45) s 

•  43L 

.  63L 

.  82L 

(0, 90) s 

i 

. 26L4 ,  . 4 3L 5 

.  63L 

(0,±45,90)s 

.47L2  ,  . 57 L  3 

.  7  7L 

.  88L 

Note:  1  No  uniform  stress  value 

2  ±45  plies 

3  0  and  90  plies 

4  90  plies 
s  0  plies 


Graphical  results  for  this  end  condition  can  be  found  in 


Appendix  A. 


Parabolic  Bending  End  Condition 

This  end  condition  reaches  a  constant  straight  line  stress 
value  for  ox  as  shown  by  the  figures  in  Appendix  B.  This 
agrees  with  the  results  Timoshenko  and  Goodier  predicted  for 
isotropic  materials.  This  constant  straight  line  stress  value 
was  utilized  to  define  the  region  of  Saint.  Venant's  principle 
validity.  As  in  the  previous  end  condition,  the  single  fiber 
orientation  plates  were  investigated  for  the  aspect  ratio  of 
five.  Again,  the  Odeg  fiber  oriented  plate  has  the  longest 
decay  length  due  to  the  relative  orientation  stiffness  prop¬ 
erties  along  the  load  axis.  The  role  played  by  the  relative 
stiffness  properties  can  be  seen  in  Fig  B-4  thru  5  for  the 
(0,  ±45,  90)s  plate.  This  is  similar  to  what  was  found  in 
the  parabolic  tension  cases.  As  in  the  previous  end  condition, 
all  of  the  decay  lengths  are  within  one  plate  width  along  the 
length  for  all  cases  considered  with  exception  to  the  (±45) s 
plate  with  aspect  ratio  of  three.  In  the  aspect  ratio  of  the 
three,  the  higher  stress  gradients  created  by  the  end  loads 
start  to  interfere  with  one  another  before  a  uniform  load  can 
be  reached.  This  caused  the  drastic  decrease  in  the  region 
of  Saint  Venant's  validity  for  the  (±45)s  plate. 

The  (0,  90) s  plate  displayed  a  much  different  behavior 
under  this  loading  than  under  the  parabolic  tension  action. 

This  is  attributed  to  the  stiffness  properties  of  each  ply. 

The  Odeg  ply  is  much  stiffer  in  tension  than  compression  due 
to  the  values  of  the  moduli,  Ej^  and  £22*  The  90ueg  ply  has 
the  same  moduli,  but  the  role  of  each  is  reversed.  Due  to 
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Poisson's  ratio  effects,  the  90deg  ply  carries  more  of  the 
load  in  this  case  than  in  the  parabolic  tension  loading. 

However,  the  stress  values  in  the  90deg  ply  is  still  less 
than  those  in  the  Odeg  ply. 

The  graphical  results  of  Appendix  B,  once  again,  depict 
variations  in  normalized  stress  values  for  different  ply 
lay-ups  due  to  differences  in  ply  stiffness  properties,  as 
all  plies  have  the  same  relative  thicknesses  for  a  given  lami¬ 
nated  plate.  The  Odeg  lamina  carried  most  of  the  stress  in 
the  plate,  as  was  found  in  all  end  conditions  investigated. 

This  is  due  to  the  increased  stiffness  properties  of  this 
orientation  in  the  direction  of  the  loads.  A  tabulation  of 
the  uniform  stress  positioning  for  this  end  condition  can 
be  found  in  Table  III. 

Clamped  End  Condition 

The  clamped  end  condition  was  investigated  to  show  con¬ 
stant  strain  triangle  modeling  difficulties  and  the  difference 
between  Saint  Venant  end  effects  and  the  moment -shear  coupling 
found  in  a  truly  clamped  end  with  rotating  clamping  devices 
and  eccentric  loading.  The  clamped  effect  is  a  displacement 
end  condition.  The  program  used  in  this  investigation  did 
not  have  provisions  for  end  displacements;  however,  the  re¬ 
sultant  force  end  condition  produced  displacements  which  agree 
with  those  reported  by  Pagano  and  Halpin.  The  clamped  end 
condition  was  examined  for  the  aspect  ratio  of  five  in  the 
single  fiber  oriented  plates.  The  laminated  plate  cases  were 
analyzed  for  aspect  ratios  of  three  and  five.  The  graphical 
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TABLE  III 

Parabolic  Bending  End  Condition 

Uniform  Stress  Station 

Orientation 

L/W  =  3 

L/W  =  5 

L/W  =10 

0 

.  76L 

-- 

90 

-- 

.  93L 

45 

-- 

.  83L 

(±45) s 

.  17L 

.761 

.  8  7L 

(0, ±45) s 

.  7  3L 

.  8  2L 

.  92L 

(0 , 90) s 

.  67L 

.  7  6L 

.  92L 

(0, ±45, 90) s 

. 56L 1  . 76L2 

.  7  6L 

.  83L 

Note:  1  ±45  plies 

2  0  aiH 

90  plies 

results  are  presented  in  Appendix  C. 

From  elementary  theory,  one  would  expect  constant  shearing 
forces  throughout  the  length  of  the  plate  and  decreasing  stresses 
due  to  moment  as  the  center  is  approached.  The  center  ox 
stress  should  be  the  constant  tension  stress  value.  The  re¬ 
sults  for  this  end  condition  generally  bears  this  out  with 
variations  due  to  material  stiffness.  The  exception  is  the 


r 


off-axis  plies.  As  displayed  in  Fig  C-3,  the  off-axis  plate 

shearing  stresses  did  not  reproduce  the  loading  condition 

curve  as  would  be  expected.  In  checking  the  average  stress 

values  at  a  given  station  along  the  x-axis,  poor  agreement 

was  found  with  the  input  load  function.  In  the  multilayered 

plate,  the  -45deg  ply  complements  the  45deg  ply  such  that, 

when  considered  together,  the  a  curve  has  a  similar  shape  as 

xy 

the  shear  end  loading  curve.  In  all  of  the  laminated  plate 
graphical  results  of  Appendix  C,  the  action  of  the  ±45  plies 
are  considered  as  acting  together  and  the  averaged  finite  ele¬ 
ment  points  are  shown  to  depict  the  stress  variations  per  ply. 

As  the  aspect  ratio  was  increased  from  three  to  five  in 
the  course  of  this  investigation,  the  variation  of  shearing 
stresses  mentioned  under  the  parabolic  tension  discussion 
increased.  The  variation  affected  the  mean  value  of  the  shear 
ing  stress  curves  from  station  to  station  as  displayed  in  the 
figures  of  Appendix  C.  The  variation  appears  to  be  related  to 
the  orientation  of  the  fibers  within  a  given  constant  strain 
triangle.  The  ±45deg  plies  were  affected  more  than  the  other 
orientations  as  the  off-axis  orientation  is  more  sensitive  to 
shear  loads.  A  possible  problem  with  modeling  off-axis  mater¬ 
ials  with  a  constant  strain  triangle  is  that  the  element  load 
conditions  are  not  the  same  as  found  in  the  continuum.  In  all 


cases,  the  a stresses  varied  with  the  shearing  stresses  as 


to  keep  the  constant  strain  triangle  in  equilibrium.  Richards 
et  al.,  experienced  variations  in  off-axis  stress  results  and 
used  averaging  over  a  very  refined  mesh  to  compensate  for  the 


variations.  Further  mesh  refinement  was  prohibited  in  this 
study  because  of  computational  storage  requirements.  The 
variation  of  the  mean  value  of  shearing  stresses  for  the 
±45deg  plies  is  on  the  order  of  ten  percent  and  for  the  Odeg 
and  90deg  plies  about  two  percent.  It  is  felt  that  this  is 
the  order  of  error  induced  by  this  phenomena.  The  ox  stresses 
did  not  exhibit  the  above  mentioned  variations.  If  variations 
were  present  in  the  c*x  results,  it  was  within  the  accuracy  of 
the  computations. 

The  practical  insight  given  is  that  the  stress  measurement 
in  off-axis  plies  can  be  drastically  different  than  would  be 
expected  from  elementary  theory,  discounting  the  variations 
mentioned  above.  This  has  implications  for  instrumenting 
laminated  plates  for  stress  measurements.  As  shown  by  the 
figures  in  Appendix  C,  the  variations  in  stress  values  per  ply 
in  a  laminated  composite  plate  can  be  considerable.  The  ex¬ 
ception  to  this  is  the  (0,  90)s  plate  shear  stress  results 
which  indicate  an  equal  division  of  the  load  between  the  plies, 
The  clamped  end  effects  are  totally  unrelated  to  the  Saint 
Venant  end  effect  and  must  be  accounted  for  in  material  pro¬ 
perty  measurements.  This  is  born  out  in  the  ±45deg  ply 
results  as  the  end  load  causes  a  constantly  changing  a 

xy 

shear  stress  in  addition  to  the  constantly  direct  stress. 

In  this  case  a  region  of  Saint  Venant's  principle  validity 
could  be  found  if  the  stress  reached  a  uniform  ^ 

distribution  and  the  o  stress  reached  a  uniform  stress 
value.  However,  no  such  station  was  found  as  the  Saint 
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stress  distributions  were  the  input 


M 


1  Venant  a  and  a 

i  x  xy 

|  functions  at  the  loaded  edge. 

(  In  all  the  figures  in  Appendices  A  through  C,  the  four 

'  stations  depicted  are  on  the  positive  half-length.  The 

percentage  of  plate  half-length  is  shown  in  each  figure. 

#  : 
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V.  Conclusions 


The  basic  conclusions  of  this  investigation  is  that 
material  properties  of  laminated  composite  plates  can  be 
measured  with  some  degree  of  confidence  in  plates  with  aspect 
ratios  greater  than  five,  but  some  striking  differences  exist 
from  one  lamina  to  another  in  a  multi-layered  plate.  The 
relative  stiffness  properties  of  a  given  lamina  greatly  affect 
the  stresses  measured.  In  some  ply  lay-up  and  end  condition 
combinations,  the  difference  in  lamina  stresses  affects  the 
decay  lengths.  These  differences  in  lamina  stresses  could 
give  rise  to  erroneous  stress  measurements  if  strain  gages 
are  employed  only  on  the  outer  lamina  of  a  plate.  In  most 
cases  with  aspect  ratios  of  five  or  greater,  the  region  of 
Saint  Venant's  validity  was  reached  at  a  distance  equal  to 
one  plate  width  from  the  edge.  Caution  should  be  exercised, 
however,  in  applying  this  result  to  other  end  conditions  as 
the  lamina  stresses  vary  greatly  with  ply  lay-up  and  end 
condition. 

In  general,  for  the  parabolic  tension  and  parabolic  bending 
condition,  the  region  of  Saint  Venant's  principle  validity  can 
be  found.  The  clamped  end  condition  points  out  the  further 
complications  incurred  in  a  tension  test  of  the  laminated 
plate.  This  end  condition  affects  laminated  plates  in  a 
manner  similar  to  the  results  reported  by  Wu  and  Thomas  for 
off-axis  plates.  The  stress  gradient  due  to  moment  does  not 
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dissipate  to  a  uniform  stress  except  at  the  center  of  the 
plate,  and  shearing  stresses  are  always  present.  If  one 
examines  the  direct  o  and  a  stresses,  a  similar  dissipa- 
tion  pattern  is  found.  The  highest  stress  gradients  are 
found  in  the  stiffen  laminas  (i.e.,  the  Odeg  ply  exhibits 
higher  direct  stresses  than  the  ±45deg  plies,  and  the  90deg 
plies  exhibit  the  smallest  stress  values) .  This  is  in  con¬ 
trast  to  the  shearing  stresses  where  the  ±45deg  plies  exhibit 
the  highest  stress  magnitudes,  and  the  Odeg  and  90deg  plies 
act  together. 

The  aspect  ratio  of  three  in  all  cases  may  not  give  a 
region  of  Saint  Venant’s  validity  depending  on  the  accuracy 
{  required.  For  this  investigation,  the  ten  percent  uniform 
stress  region  was  only  a  small  region  in  the  center,  and  for 
uniformity  of  stress  less  than  ten  percent,  there  is  no 
region  of  Saint  Venant's  validity. 

More  can  be  drawn  from  this  study  if  one  thinks  of  the 
possibility  of  superimposing  the  stress  fields  produced  by 
the  three  stress  fields.  This  would  model  a  non-uniformly 
clamped  end  with  a  rotating  grip.  In  this  case,  the  Saint 
Venant  end  effects  would  decay  at  about  one  plate  width  from 
the  loaded  edge,  but  the  stress  field  due  to  the  moment-shear 
coupling  effects  would  not  dissipate.  The  simple  load  per 
unit  area,  P/A,  stress  field  would  only  be  realized  at  the 
center  of  the  plate,  regardless  of  aspect  ratio.  In  addition 
to  the  P/A  stress  field  for  ox,  a  parabolic  shearing  stress 
would  be  present.  This  leads  to  the  conclusion  that 
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meaningful  ax  stress  measurement  can  only  be  made  at  the 
center  of  the  plate  in  material  property  calculations. 
Shearing  stress  measurements  would  be  greatly  influenced  by 
the  shear-moment  coupling.  In  addition,  the  stress  gradient 
per  ply  orientation  are  not  the  same.  This  can  lead  to  erro 
neous  stress  measurements  if  only  the  outer  plies  are 
instrumented . 

It  is  further  concluded  that  not  enough  data  has  been 
gathered  to  adequately  investigate  the  variation  in  stress 
values  found  in  the  off-axis  plate  cases  and  in  the  clamped 
end  condition  where  shear  loading  is  used.  The  effects 
of  using  constant  strain  triangles  with  off-axis  orthotropic 
materials  needs  further  study.  It  is  felt  that  further  mesh 
refinement  and  averaging  stresses  across  several  triangles 
would  produce  better  results;  however,  computer  storage 
limitations  are  a  problem  with  the  program  utilized  herein. 
This  aspect  of  finite  element  analysis  merits  further  study. 


I 


I 

I 

i 


Bibliography 


1.  Timoshenko,  S.P.  and  J.N.  Goodier.  Theory  of  Elasticity. 

Third  Edition.  New  York:  McGraw-Hill  Book  Company, 

1970. 

2.  Ashton,  J.E.,  J.C.  Halpin,  and  P.H.  Petit.  Primer  on 

Composite  Materials:  Analysis.  Westport,  CT : 

Tecnnomic  Publishing  Co.,  Inc .  ,  1969. 

3.  Popov,  E.P.  Introduction  to  Mechanics  of  Solids. 

Englewood  Cliffs :  Prentice-Hall ,  Inc .  ,  1 963 . 

4.  Horgan,  C.O.  "The  axisymmetric  End  Problem  for  Transversely 

Isotropic  Circular  Cylinders."  International  Journal  of 
Solid  Structures,  10:837-852  (1974)  . 

5.  Folkes,  M.J.  and  R.G.C.  Arridge.  "The  Measurement  of 

Shear  Modulus  in  Highly  Anisotropic  Materials:  The 
Validity  of.  St.  Venant's  Principle."  Journal  of 
Physics  D:  Applied  Physics,  8^:1053-64  (1975). 

6.  Choi,  I.  and  C.O.  Horgan.  "Saint  Venant's  Principle  and 

End  Effects  in  Anisotropic  Elasticity."  Transaction 
of  the  American  Society  of  Mechanical  Engineers: 

Journal  of  Applied  Mechanics:  424-4~30  (September  1977)  . 

7.  Pagano,  N.J.  and  J.C.  Halpin.  "Influence  of  End  Constraint 

in  the  Testing  of  Anisotropic  Bodies."  Journal  of 
Composite  Material,  2_:  18-31  (January  1968)  . 

8.  Jones,  R.M.  Mechanics  of  Composite  Materials.  Washington, 

D.C:  Scripta  Book  Company,  1975. 

9.  Whitney,  J.M.  Testing  and  Characterization  of  Composite 

Materials.  Technical  Report  AFML-TR- 71-124 .  Wright  - 
Patterson  AFB,  OH:  Air  Force  Materials  Laboratory, 

June  1971. 

10.  Wu,  E.M.  and  R.L.  Thomas.  "Off-Axis  Test  of  a  Composite." 

Journal  of  Composite  Materials,  2_:  523-56  (October  1968). 

11.  Rizzo,  R.R.  "More  on  the  Influence  of  End  Constraints  on 

Off-Axis  Tensile  Tests."  Journal  of  Composite  Materials, 
3:202-19  (April  1969)  . 

12.  Richards,  G.L.,  T.P.  Airhart,  and  J.E.  Ashton.  "Off- 

Axis  Tensile  Coupon  Testing."  Journal  of  Composite 
Materials ,  3:  586-89  (July  1969). 


37 


13.  Fung,  Y.C.  Foundations  of  Solid  Mechanics.  Englewood 

Cliffs:  Prentice-Hall ,  Inc.,  1965 

14.  Filonenko-Boradich,  M.  Theory  of  Elasticity.  New  York 

Dover  Publications,  Inc .  ,  1965 . 

15.  Sokolnikoff,  I.S.  Mathematical  Theory  of  Elasticity. 

New  York:  McGraw-Hill  Book  Company,  1956 . 

16.  Pipes,  R.B.  and  N.J.  Pagano.  "Interlaminar  Stress  in 

Composite  Laminates  Under  Uniform  Axial  Extension." 
Journal  of  Composite  Materials,  4:538-548  (October 
1970) . 

17.  Durocher,  L.L.  and  A.N.  Palazotto.  "Notes  For  Computer 

Methods  in  Stress  and  Structural  Analysis."  Lecture 
material  used  in  MC6.41,  Fundamentals  of  Matrix 
Structural  Analysis.  School  of  Engineering,  Air 
Force  Institute  of  Technology,  Wright-Patterson  AFB, 
Ohio . 

18.  Venkayya,  V.B.  Finite  element  program  OPSTAT,  written 

at  the  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Patterson  AFB,  Ohio. 


38 


Appendix  A 


Graphical  Results  of  the  Parabolic 
Tension  End  Condition 

The  following  figures  resulted  from  the  use  of  the  quarter 
plate  mesh  shown  in  Fig  4-C  with  the  exception  of  Fig  A-3, 
which  utilized  the  full  plate  mesh  of  Fig  5-b.  All  stresses 
in  the  figures  are  normalized  to  the  average  stress  oQ  as 
follows : 

a 

SIGMAX  =  (A-l) 

where  SIGMAX  is  the  normalized  stress  value,  ox  is  the  stress 

in  the  x~direction  and  aQ  is  the  average  stress  computed  from 

N  the  average  resultant  force.  In  all  figures  the  four 
o 

stations  depicted  are  on  the  positive  half  length  of  the  plate. 
The  percentage  of  the  half  length  is  noted  in  the  legend  for 
the  station  plotted. 
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Appendix  B 

Graphical  Results  of  the  Parabolic 
Bending  End  Condition 

The  following  figures  resulted  from  the  use  of  the  full 
plate  mesh  shown  in  Fig  6-A.  All  stresses  in  the  figures 
are  normalized  to  the  average  stress  oQ  given  by  Eq  (A-l). 
Legend  comments  found  in  Appendix  A  apply. 
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Appendix  C 


Graphical  Results  of  the 
Clamped  End  Condition 


The  following  figures  resulted  from  the  use  of  the  full 
plate  mesh  shown  in  Fib  6-B.  All  stress  in  the  figures  are 
normalized  to  the  average  stress  oq  given  by 

a 

SIGMAX  =  — 
o 


and 


SIGMAXY  = 


(C- 1) 


where  ax  is  the  stress  in  the  x-direction,  axy  is  the  shearing 
stress,  SIGMAX  is  the  normalized  ax  stress,  SIGMAXY  is  the 
normalized  axy  stress,  and  aQ  is  the  average  stress  computed 
from  Nq ,  the  average  resultant  force.  Legend  comments  of 
Appendix  A  apply. 
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Analytic  Solution  for  the 
Orthotropic  Plate 

The  following  mathematical  development  follows  the  analytic 
solution  presented  by  Choi  and  Horgan.  An  approach  similar  to 
this  can  be  followed  if  laminate  plate  material  properties  are 
considered  in  place  of  constitutive  constants.  The  development 
for  the  Odeg  ply  is  presented  here  for  completeness.  For 
further  details,  Reference  6  should  be  consulted. 

The  strains  for  a  plate  in  a  state  of  plane  strain  are 
as  follows: 


2e. -  =  U.  .  +  U.  . 

13  i,3  3,i 


(D-l) 


where  y  is  the  displacement.  The  strains  can  also  be  written 
in  terms  of  the  in-plane  stresses  as  follows 


ell  = 

enT 

11 + 

812t22  + 

B16t12 

e22  " 

B21t 

11 + 

B22t22  + 

626t12 

2e12  * 

661t 

11  + 

B62T  22  + 

B66t12 

(D-2) 


where  the  subscripts  on  the  stresses  and  strains  represent  the 
planes  of  symmetry  in  the  material  and  the  3’s  form  the  com¬ 
pliance  matrix.  The  elastic  constants  (g’s)  are  assumed  to  be 
such  that  the  strain  energy  density  is  positive-definite,  and 


The  Airy  stress  functions  used  are 

a=(b  .  a  *  <b  and  t  *  -d>  fD-4) 

x  ¥yy  *  y  ^xx  xy  Yxy  v  J 

where  the  subscripts  denote  differentiation. 

The  compatibility  relations  yield 

^22^xxxx  ‘  2B26*xxxy  +  C2612  +  666)4,xxyy 

'  2016*xyyy  +  ^ll^yyyy  =  0  (D*5) 

The  equilibrium  conditions  are  identically  satisfied  by 
the  stress  functions.  A  solution  of  Eq  (D-5)  is  sought  in 
the  form 

*  -  e*YXF(y)  (D-6) 

where  y  is  a  constant. 

The  lateral  boundary  conditions  are 

F(±l)  =0  F* (±1)  =  0  (D-7) 

for  a  plate  loaded  in  tension. 

Substituting  the  above  stress  function  into  the  compati¬ 
bility  relation  (D-5)  yields: 
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^llF,yyyy  +  2^16YF,yyy  +  (2ei2  +  e66^Y  F, 


yy 


*  «26^F,y  +  822''4p  =  0 


CD-8) 


where  the  subscripts  on  F  indicate  differentiation. 

Eqs  (D-7)  and  (D-8)  form  an  eigenvalue  problem  for  the 
ordinary  differential  Eq  (D-8).  A  solution  of  Eq  (D-8)  is 
sought  in  the  form  of 


or 


F(y)  =  ewy 


<J> 


e-YXe“y 


=  e.wy'YX 


(D-9) 


(D- 10) 


Eq  (D-8)  becomes 

8nu>4  +  2816yo)'5  +  (2B12  +  866)y2o)‘'  +  2B26y3w 

♦  B22Y4  =  0  CD-11) 

We  now  let  to  =  yY  where  to  and  y  are  complex  numbers. 

Eq  (D- 11)  becomes 

glly4  +  2g16y3  +  ^2012  +  ®66^y2 

+  2B26y  +  e22  =  0  (D’12) 


From  the  positive-definite  strain  energy  density  conditions, 
it  can  be  shown  that  the  roots  of  y  are  complex  or  purely 
imaginary. 


The  conjugate  pairs  of  roots  are 


“1,2  =  Cpl  4  ^1)Y  63,4  =  ^P2  4  q2)Y  CD'13) 


where  p  =  P  ±  iq  (D-14) 

Using  the  Eqs  (D-9)  and  D-13),  the  complete  solution  is 


(Pi±iqi)Yy  (P^  +  iqOyy 

F(y)  -  Axe  1  4  +  A2e  L  1 


CD-15) 


PiYy 


1  l  T  / 

F(y)  =  e  (C1  cos  q^y  +  C2  sin  qjyy) 


P2  y 

+  e  (C3  cos  q2  y  +  Cj  sin  q3  y) 


CD-16) 


where  and  are  arbitrary  constants. 

For  the  orthotropic  case,  B16  =  0  and  B26  =  0  and  three 
cases  of  the  solution  CD -16)  can  be  considered. 


Case  A: 

(2^2  + 

* 

4611B22  > 

Case  B: 

C2Bi2  ^ 

he)2  - 

4B11B22  < 

1  Case  C: 

(26l2  + 

he)2  - 

4eilS22  = 

For  Boron  Epoxy: 

el2 

*  .227 

i 

i 

he 

*  24.2 

hi 

=  .722 

622 

*  13.6 

(D- 17) 


(26i2  +  866)  -  46116 


1122 


=  [2  C -  227)  +  24 . 2  j  -  4  (.  722) (13 . 6)  >  0 


•*.  Case  A  applies 


when  Case  A  holds 


~  0  ,  w  is  purely  imaginary 


“1,2  =  ±qlY  ’  “3,4  =  ±q2Y 


The  even  eigenfunctions  of  Eq  (D-8)  are 

Fe(y)  =  cos  -  cos  Yq2Y 

cos  yq^  cos  yq2 

where  y  is  a  root  of 

qx  tan  yqx  -  q2  tan  Yq2  =  0 

The  odd  eigenfunctions  of  Eq  (D-8)  are 


Fe  (y) 


where  y  is  a  root  of 


sin  yq xy  sin  yq 2y 

sin  yq.^  sin  yq2 


(D- 18) 


(D- 19) 


q2  cot  yq2  -  q2  cot  yq2  =  0 


From  the  Airy  stress  functions  and  Eq  (D-6) ,  we  get 


°x  =  e"YXF"  °y  =  Y2e'YXF  x  =  +ye  *yxF 1 


(D- 20) 


The  total  stress  can  be  found  by  taking  a  summation  of 


the  individual  solutions  or 


$  =  Z  C,e'YkxFv(y) 


^  '-V 

k=l  k 


aY  =  Z  Cke*YkXFk”  (y) 


k=l 


(D- 21) 


°xy  "  ’k‘ Fk^ 


Txy  ■  ./Ar^k^Cy) 


k=l 


where  Ck  is  an  arbitrary  constant  determined  from  the  end 
condition.  Using  the  relationship 

ez  =  2(cosh  z  +  sinh  z)  (D-22) 


Eq  (17)  becomes 

* 


Z  (ak  cosh  Yk  +  bk  sinh  Ykx)Fk(y) 
k  1 


a  = 


a..  = 


E  (ak  cosh  Ykx  +  bk  sinh  Ykx)F£(y) 


I  Yk(ak  cosh  Ykx  +  bk  sinh  Ykx)Fk(y) 
k  1 


(D-23) 


xy 


Z  y(ak  sinh  Ykx  +  bk  cosh  Yj.x)F£(y) 
k  1 


where  ak  and  bk  are  arbitrary  constants  determined  from  the 
end  conditions  and  orthogonality  conditions. 
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